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ON ■ Abstract 
^ . 

We analyze new QCD sum rules for the pion wave function (WF) ipj^{x), obtained recently in 
I in the non-local vacuum condensate method for non-diagonal correlators, and suggest a new 

approach for extracting WF of the 7r-meson and the mass and WF of the first resonance. As a 
result, we obtain approximately the same form of the pion WF by two different methods and for 
two different ansatzes of non-local quark condensates. We predict the mass of the 7r'-resonance 
' and obtain the form of its WF. 

> 

m ■ 1 Introduction 

^ I An important problem in the theory of strong interactions is to calculate, from the first prin- 
OS I ciples of QCD, the distribution functions fp/nix) and hadronic wave functions . . . 

ipN{xi,X2,X3), ... which accumulate all the necessary information about non-perturbative 
long-distance dynamics of the theory. These phenomenological functions appear naturally as 
Q-i! a result of applying "factorization theorems" to hard inclusive and exclusive processes [|], 0], 

Another kind of such phenomenological numbers are the quark and gluon condensates (: 
g(0)g(0) :), (: G(0)G(0) :)..., the basic parameters of the QCD SR approach reflecting 
the non-perturbative nature of the QCD vacuum. A usual practice is to calculate the hadronic 
functions /(x), by using the condensates as input parameters p. 

The situation with hadronic WF's is more complicated than with distribution functions: 
the firsts appear only in an integrated form like some convolutions. It seems now that only the 
QCD Sum Rule (SR) approach and lattice calculation [^] can provide an information about 
the form of hadronic WF. The most popular set of hadronic WF, due to V. L. Chernyak, 
A. R. Zhitnitsky and I. R. Zhitnitsky (CZ) p|, was produced with the help of QCD SR for 
the first moments of WF's. These SR were based on the diagonal correlator of appropriate 
axial currents and the condensates of the lowest dimensions; the pion WF thus obtained has 
the well-known "two-hump" form. But now it is known that the hadronic functions are rather 
sensitive to the structure of non-perturbative vacuum Therefore one should use a non-local 
condensate like (: q{0)E{0, z)q{z) :) which can reflect the complicated structure of the QCD 
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vacuum. (Here E{0,z) = P exp{i Jq^ (it^y4jj(t)ra) is the Schwinger phase factor required for 
gauge invariance.) 

Earher, one of the authors (S. M. together with A. Radyushkin) constructed a modified SR 
with non-local condensates and demonstrated that the introduction of the correlation length 
A for condensate distributions produces much smaller values for the first moments of pion WF 



than the CZ values [|T0|, [|TT1- This leads to the form of the pion WF strongly different from the 
CZ form and close to the asymptotic form if'^^{x) = 6x{l —x). Our goal in the present paper is 
to obtain directly the form of the pion WF and the first resonance, using the available ansatz 
for non-local condensates (section 2). This program has been suggested recently and realized 
in [|^, we develop alternative methods of extracting wave functions from this sum rule for a 
non-diagonal correlator. We suggest two different ways: the first uses integrated (in a sense) 
properties of ansatzes (see section ^) and may be considered as an alternative to traditional 
Borel SR; while the second uses their local properties (see section ^ and follows to the Borel 
SR idea. The results obtained by these independent methods agree with each other. 



2 Non-diagonal correlators and pion WF 
2.1 The models of non-local quark condensates 

It is appropriate to recall some general features of derivation of the QCD sum rules with non- 
local condensates 0, 0], [|l6l. The quark and gluon fields are taken in the Fock-Schwinger 
gauge A^{z)z^ = 0, therefore path-ordered exponentials E{0,z) = Pexp{i /q^ (it^y4^(t)ra) are 
equal to 1. It is convenient to parametrize the z^-dependence of the simplest bilocal quark 
condensate (: q{0)q{z) :) = (: g(0)g(0) ■.)Q{z'^) in analogy with the a-representation of a 
propagator R 

Qiz')= e^^'l'f{s)ds. (1) 







The correlation function f{s) may be interpreted as "the distribution function of quarks in the 



vacuum" since its nth moment is proportional to the matrix element of the local operator 
with to nth power: 

-.vw<;.^^^f^. (2) 

r(n + 2) (:gg:) 

For the lowest two moments one obtains 

f(\^ 1 r f(\^ 1 (: qD\ :) _ 

f{s)ds = l] / sf{s)ds = - = — , (3) 

Jo 2 {: qq :) 2 



with Aq meaning the average virtuality of vacuum quarks. 



2 
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In a similar way, parametrizing the quark-gluon non-local condensate 

(: mwicrGiOMz) :) ^ (: g(0)g(0) :)Q,{z'), 



^ In deriving these sum rules we can always make a Wick rotation, i.e., we assume that all coordinates are 
Euclidean, < 0. 
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one can introduce the quark-gluon distribution function fi{s). (Note the approximate character 
of this definition: its l.h.s. appears as a three-point correlator and only after reducing to two- 
point geometry can be parametrized in terms of two-quark non-local condensates [|10|] .) 

To construct models of non-local condensates, one should satisfy some constraints. For 
instance, if we assume vacuum matrix elements (: q{D'^Yq :) to exist, then the function f{s) 
should decay faster than as s ^ oo. If all such matrix elements exist (for all /c), a 

possible choice could be a function /(s) ~ I" , or /(s) ~ e'^l" ^ etc. at large s. The 
opposite, small-s limit of /(s) is determined by the large-|2| properties of the function Q(-2^). 
This behavior has been analyzed in detail in in the framework of QCD SR for heavy quark 



effective theory (HQET): it has been demonstrated that for a large Euclidean 2;, non-local 
quark condensate Q(z^) ~ e~'^'^ with A = [Mq — mQ)|mQ^oo being the lowest energy level of 
the mesons in HQET (numerically, A is around 0.45 GeV). This means that /(s) ~ e^^^^* 
in the small-s region. Another hint to prove this behavior of f{s) can be obtained from the 



results of the lattice calculation [T^ where the exponential decay behavior for the correlator 
(: G{0)E{0, x)G{x) :) of the gluon vacuum strengths has been found. For this reason one should 
expect in a theory with confinement 

(:g(0)g(z):) ~ exp(-|z|A) (4) 

where A is the correlation length. This law generates the above-mentioned behavior for f{s) 
at small s. So, we arrive at the set of ansatzes AI: 

Al : f{s) = Ni ■ exp A^/s — s'^af^ — Gaussian decay 

A2 : f{s) = A^2 ■ exp (— A^/s — S(J2) — exponential decay (5) 



AO exp(-A2/s)l /AlV'"'^ 

A3fe : fk{s) = Y{k-1) s VV i ^ inverse power decay 

where A^ = 0.2 GeV"^ and A^ = A^(/c — 2), the normalization constants iVj and the cxj-parameters 
are fixed by eq.(^), where for the average virtuality of vacuum quarks we take the usual QCD 
sum rule value Il| ~ OAGeV^. 

Ansatz Al has been used in it is effectively close to ansatz A2. The second ansatz leads 
to the same behavior for the (: q{0)q{z) :) as the massive causal scalar propagator with a shifted 
argument does: 

/yrA 

(^2 + \^^y 



{: mq{^) -K - ^2-— ^r^^— ^exp(-2AVa2+k2|/4) (6) 



-00 



A^2 = — I 7 ^li^) is the modified Bessel function. 



2^K, (2AV^ 



The ansatz A3 has analytic properties in s different from those of A2, is very convenient for 
calculations and imitates the 5(s— A) when the parameter A; ^ 1. But it doesn't look sufficiently 
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"physical" due to the absence of higher local condensates mentioned above. Nevertheless, we 
shall test it also at the values k = 3, 4, 5. This ansatz leads to the asymptotic behavior 



2fc-3 



: mii^) ■■) . - T^TT^ ('4) ' exp {-A\z\) . (7) 



r(A;- 1) V 2 
2.2 Sum rule 

Let us write down the model sum rule for the "axial" wave functions y^^... (x) of the pseudoscalar 
mesons suggested in [|I|: 

+ ^^,{x)e-'^l'"'' + ^A^)e-<""'' + ...^^{^^,x) (8) 

with the function /(s) specified in the preceding section. It should be noted that this sum rule 
results from the approximations (for a detailed discussion, see [jl[): 

1. reduction of three-point correlators to two-point ones; 

2. /i(.) = 2A2./(.); 

3. contributions of higher non-local condensates {qGGq), {qGGGq), . . . could be neglected. 

As has been shown in the function $( , a;), i.e., the weighted sum of all WF's, is given 
by two humps (one centered at xa_ = sa/M"^, where sa is the point of maximum for the ansatz 
distribution function /(s), and the other - at Xyi = 1 — sa/M"^) moving as changes. When 
increases, the humps become narrower, higher and more close to boundary points x = or 
X = 1. For = 1 GeV^, e.g., the function $(-^,x) looks very much like the CZ WF (see 
Fig.|l]). However, $(jg5-,x) is not just the pion WF: the larger M^, the larger contamination 
from higher states. Moreover, the ("axial") WF of the higher pseudoscalar mesons ir', ir", . . . 
produce zero total integrals (whereas the pion WF is normalized to unity, see also |]T^) and 
therefore they should oscillate (see, e.g., Fig.^) 

At low M^, the pion WF dominates in the total sum ^{j^,x) (however, one cannot take 
too low because the operator product expansion fails for < A^). When = 0.4 GeV^, 
it was observed that $(-^,a;) is very close to the asymptotic wave function of the pion (see 
Fig-D- 

In our opinion, this very nice picture of successive switching on resonances has an essential 
drawback: there is no strict criterion for selecting the value of for determining WF of the 
ground state (pion). The author of ^ writes: "... one should be more accurate here, since 
even the modest increases of to 0.5 GeV^ or 0.6 GeV^ induce humps in $(x,M^)" (see 
Fig-D . 
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3 Method of Integral Projector 



As we have seen before, all the resonances contribute to ^{-^,x) and the problem is to dis- 
tinguish between the ground state {f{x)) and resonances {(pi{x)) and between resonances. We 
should know the whole spectrum /ij to extract the ground state. We represent here a method 
to obtain ipi{x) and fii step by step. It is based on an integral transformation of equation (H). 

As a result, we shall obtain approximately the same form for the ground state and the first 
resonance by using essentially different ansatzes A2 and A3. 



3.1 General formulas 

In what follows it's convenient to use a new variable r = instead of the Borel parameter 

M^. The r.h.s. of SR (W), i.e. $(r, x) is defined for r G [0, yj]- Let us consider this function in 

the whole complex plane of r and define the projector P{N, tq, uj): 

where the contour C is a vertical line (7 = (c — zoo, c + zoo) with c > l/A^, the c lyes righter 
than any pole of the integrand, u > and tq G (0, ya]. To obtain the result of this projector 
action on the initial SR-representation 

V^i^) + exp(-r/ii) = $(r, x), (10) 

1=1 

consider its action on a simple exponential. Evidently, by closing the contour C to the left we 
get (due to the well-known residues theorem): 

P{N,To,uj)exp{-^T) = 6{uj - fi) ■ exp{-fiTo) il - -j . (11) 



Then, for > 0, we have: 



UJ 



(p{x) - ('^ -) exp{-Tofi,) = P{N,to,uj)^{t,x) = ^n{oo,x). (12) 



To understand this result, it is instructive to consider the case N = 0, tq = 0: 

ip{x) + ^ii^) = ^(0, 0, u;)$(r, x) (13) 

i>l 

We see that by varying u one can switch on (switch off) more and more resonances in the 
l.h.s. of SR. For > one can successively determine the positions (masses /ij) of resonances. 
Formula (|T^ ) is the main result of this subsection. 

Consider now the meaning of the introduced parameters of the projector P{N, Tq, u): 
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1. The parameter uj has a clear physical meaning — it produces a real division of the set 
of resonances ( with fii > uj). Contrary to the role of the "unphysical" Borel parameter 

that supplies an exponential cut of higher resonances, it gives a (smoothed by powers 

( 1 j ) real step cut. 

2. The real life is more complicated than the simple model (plsl) : to suppress the non-desirable 
asymptotic behavior of the function $(r, x) for r — >■ oo (in order to neglect the integral 
over a semicircle closing the initial contour C), one should use the parameter > 0. 
Besides, the projected SR with different can be used for determining the subsequent 
resonances. 

3. At last. To defines the significant region of variable r which mainly contributes to the 
integral. 

Let us consider the third feature in detail. Note that the physical meaning has only large- 
\z\ asymptotics of the non-local quark condensate (: q{0)q{z) :). But no information about 
this behavior can be extracted from the known values of the first few moments. Contrary to 
that, the behavior of the condensate at small \z\ isn't fixed at all, and different kinds of the 
behavior of the correlation function f{s) are allowed at a large s (or a small r). The reason 
is the dominance of perturbative contribution in this region. Different ansatzes Al, A2, A3 
have been suggested earlier in this way. It is easy to show that (^ generates the exponential 
behavior of the correlation function f{s) at a small s, f{s) ~ exp(— A^/s) that corresponds to 
the law 

$(r, x) ~ exp(-rA^) 

at a large r for the r.h.s. of SR. 

As the behavior of $(t) at small r isn't fixed we always choose the value of tq in some 
region near point Since the final result, in any way, should not depend on the behavior 

of $(t) at a small r <^ this is a constraint on the available correlation function 

/(1/r) and the value of tq. The practical recipe is that one takes the parameter tq so that 
the contribution of the region near point r = to the integral (P) must be small (see the 
integration with ansatz A2 in the next subsection). 

Consider now as an illustration, the result of the action of P{N, tq, w) on the n-th term of 
the operator product expansion (OPE) (0„) -r" in the r.h.s. of the traditional SR, proportional 
to r". After explicit integration we obtain for N > n: 

One can see the re-formation of the usual set of power corrections: 

In the limiting case ujtq ^ oo (at a fixed Tq), when all the resonances are included into the 
l.h.s., the projector P{N,to,u) is similar to the projector on a single value of r, i.e. its kernel 
in the integral representation is close to 5(r — Tq): 
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and the SR is restored again to its initial form (with a new argument r —>■ tq). 

For intermediate values of uj, the cut of the tail of resonances in the l.h.s. (i. e. those 
which satisfy uj < fii) leads to the relative growth of theoretical contributions in the r.h.s. of 
the SR (see (0)) with respect to the remaining resonance contributions. In other words, it can 
be interpreted as an effective increase in the scale of the standard vacuum expectation values 
(On) ■ This example demonstrates, in particular, the strong dependence of SR treatment on the 
adopted model for the phenomenological part of SR. 



3.2 How to extract a ground state WF 

Let us assume that the experimental mass values are known /io = 0, fii = 1.7, ... and they 
correspond to ansatzes A2 or A3. Then, we can obtain an expression for v^(x) by employing 
equation ([T^ ) with the parameter uj = ooi < fii that corresponds to the saturation of the ground 
state (p{x) 

ip{x) = P{N, To, tui)$(r, x) = ^Niooux). (15) 

Then, the final result should not depend on the parameter strongly. However, in the real 
life we don't know the position fii exactly, therefore we may get a contamination from the next 
state. This contamination becomes smaller when A^ increases, due to the power suppression 
near threshold /ij. So, the following criterion of saturation of the ground state can be suggested: 
the form of the curve ip{x) = doesn't depend on A^ at some Ui. This is a very simple 

but no so effective criterion for the following reason: the region of comparatively small u (i. 
e. uj < 1 GeV"^) seems to be transient, here the corresponding WF is not well normalized 
(its normalization is of the order of 0.5 as compared to needed 1). The normalization of WF 
restores at > 1.5 GeV^ , but the correction for the first resonance is too small to be extracted 
with the help of this criterion. 



To integrate the r.h.s. of SR ([15]), we take residues: there are two residues in the case 
of ansatz A2 at the points t = tq and r = (essential singularity). At the end we obtain 
that the saturation realizes near u; = 1.8 -r- 2.0 GeV"^ and the function $Ar(u;, x) doesn't change 
noticeably for A^ = 0, 1, 2, 3. The parameter tq is chosen in accord with the suggested criterion 
of smallness of the contribution from the residue at r = and appears to be tq ~ 2.5 GeV~^. 
The curve is in form rather close to the asymptotical WF, and the second moment of WF 
(^^) =0.19 (the corresponding curves are shown in Fig.|^). 

We have completely different behavior of the r.h.s. of SR near r ~ for ansatz A3fc, 
$(r, x) . For this reason, we are should use the projector with the parameter N = k + 2 
for this ansatz. We proceed here with ansatzes with k = 3,4,5. The integral in the r.h.s., 
i>7v(<^0!^)) is determined by the one residue at the point r = tq. Performing the procedure 
like for ansatz A2 with the same value of Tq, we obtain slightly thinner forms of (p{x) for all 
these ansatzes. As it will be clear later, the most adequate to pion physics is ansatze ^3^ with 
3 < k < 4. For this case we observe the saturation of the r.h.s. of SR at = cuq ~ 2 GeV^ 
with (^^) ^ 0.16. (We could conclude about this saturation by the smallness of variations of 
^n{^^o,x) with A^, of course the WF normalization now reaches unity.) The results are also 
close to the asymptotic WF, see Fig.|3. 
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3.3 How to extract resonance WF's 



Suppose we know the values of /ii, /i2, . . ., then we can extract the first resonance by inverting 
equation ([12|) , e.g., at iV = 1: 



('li(^,x)-l>i(^i,x) 

^1 w = ^^^^^^ ^ ^ 

where /^i < < /i2 corresponds to saturation of the first resonance. 

The later equation may be improved if one takes into account a more detailed spectrum 
model. Let us suggest that the resonance contribution to the spectral density has a finite width 
2A, e.g.: 

6{s — jj) ■ ^i{x) — s> p(s — yU, A) ■ where p{s — /i, A) |a-»o = ^(-^ — jj) 

This leads to a more complicated expression in the l.h.s. of SR than simple exponentials. 
Applying the projector P{N, tq, uj) to the l.h.s. of SR with this new spectral density p one can 
easily obtain 

poo I d \ 1'°^ 

P(N,To,Lj) e-'^p(s- p,A)ds= [1 + ^—] / e(u-s)e~'^''p(s-p,A)ds (17) 

JO \ udTQ J JO 

instead of expression (|n|). To estimate the main effect of the finite width A, we apply the 
simplest (and rather rough) "step-model" for p 

p{s -p,A)=9{fi + A> s)e{s >p-A)-^, (18) 

at natural constraint A <^ p. Substituting (|TBp into (|r^) and integrating we obtain a new 
expression for the l.h.s. of SR. The latter leads to modified equation for ipi{x)] at iV = 1 and 
> /i + A we arrive at the expression 



$i(cj,2;) - $i(/ii,x) Arn / 1 \ 

= ^ -P (-0.0 ■ A,, = A . (coth(A.„) - — ) (19) 

The final numerical result due to finite width in ([T^) is not large for parameters of the first 
15| A ^ 0.18 GeV^, pi ~ l.lGeV'^ but may be important for the next resonance. 



resonance 



Let us assume now that we don't know the spectrum {pi}. How can we find a resonance for 
the particular ansatz? Evidently, equation ([T^ has the best sensitivity to it at the lowest A^, 
because the threshold effects are suppressed for large N . However SR with = is not also 
good: there is a too rough simulation of the continuous (over uj) r.h.s. by the set of step-like 
contributions. 

The case of SR with A^ = 1 is more attractive: here steps are smoothed by switching- 
functions {l — and the discontinuity is present only in derivatives of the l.h.s. of SR. As 
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it has been mentioned in the discussion of the criterion for selecting resonances in the former 
subsection, the magnitude of this discontinuity is too small to be trigger resonances. If we 
examine the difference x) — '^^^1(00, x), then we obtain the very convenient criterion for 

resonances: 

9 {cj - Hi) (fi{x) ■ exp{-niTo) — (1 ^1 = ^n{(^, x) - ^N+iiuj, x). (20) 

UJ \ UJ / 

As one can see, at = u;* = /^i, this difference is equal to zero and for uj > [i\ becomes negative, 
if f^iix) < 0, otherwise, positive, with l/cu-decreasing at large u. For ansatz A2 we see in Fig.H 
that the zero of this difference atA^=l, x = 0.5is reached for a;* = 1.8 GeV'^ that is quite 
reasonable for the pion case (the experimental value is fiiexp ~ 1-7 ± 0.25 GeV'^). It should be 
emphasized that $1(0;, x) imitates the l.h.s. of SR as a function of u rather well despite of their 
different origin. Note here that the position of the root uj^ depends also on the value of x: in 
the region 0.35 < x < 0.65 this dependence is rather weak. 

So, we have found the position of the first resonance fii that agrees with ^lexp, based 
on ansatz A2. This is an unexpectedly good result for our rather crude model of non-local 
condensates (see subsection |2.2|) . The saturation by the first resonance is reached near uj ^ 
2.5 ^ 3 GeV^. The curves corresponding to the resonance (pi{x) (for different values of uj near 
UJ = 3 GeV"^) are shown on Fig.|. 

The same program for ansatz A3 with k = 3,4,5 produces the roots uj^,{k) = 1.3, 3.0, 
5.5 GeV^, respectively. For this reason the choice 3 < < 4 is more adequate to the pion 
case. Here we also observe the saturation of normalizations and forms of the ground state WF 
before switching on the first resonance (i.e. at uj < uj^,{k)). The curves of the ground state WF 
corresponding different A; = 3, 4, 5 are very close. 

4 Method of Differential Projectors 
4.1 General formulae 

We shall construct here a trick to extract the ground state and resonances, based on the 
criterion traditional for QCD SR approach: a physical quantity, (p{x), should not depend on 
the Borel parameter on the stability plateau. However, as we have seen in section 2, the 
weighted sum of all WF in the initial representation (^ has a strong dependence on even 
for < 1 GeV^ due to contamination by the resonances, and there is no stability of the form 
of a curve at all. 

Let us suggest that the positions of the first n resonances are known exactly. Then, one 
may hope to obtain the ground state ip{x) that doesn't depend on in the region [A^,M^], 
by cancelling the contributions of these resonances. The contamination in this case will be 
defined by the contributions of the (n + l)-th resonance and subsequent resonances and will be 
suppressed by the exponential factor ~ exp{—fin_^_^/M'^). If the positions of some resonances 

^ It is similar to the formulation of the "inverse problem": to obtain the (p{x) one uses the known spectrum 
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are not well known, one can vary the parameters /i, to obtain the best fit of stability. Of 
course, this proposal will "work" if the resonances are well distinguished and the corresponding 
pre-exponential factor ipn+i{x) is not too large (for any x). Note at least that the renorm-group 
evolution of the condensates with respect to should be taken into account to complete this 
analysis in a separate publication. 

To cancel the contribution of a set of resonances, we shall use a differential operator. Ap- 
parently, one can obtain for one resonance (fi ( here D = dr )'■ 

(^l + e~^^^ipi{x) = 0; (^1 + ^(x) = ifix). 

From this result it is easy to find the action of the projector Qo{n) f\ 

Qo{n)^fl(l + -], (21) 



i=l 



on a finite set of the resonances L„ = exp(— /ijr)v9j(x): 

Qo{n) fEe-^'>.(x)^ =0; Qo{n) ( ^{x) + j2e-^'^V^{x)] = ^(x). (22) 



^i=l / \ i=l 



So, applying the projector Qo{n) to eq. (p!0D, one obtains the expression for the ground state 
ip{x) 

ifix) = QoinMr, x) + O(e-'^"+^0- (23) 

We can extend the region of stability for ip{x) with respect to r (or M^) by adding subsequent 
resonances to projector Qo{n) (and making the last term in eq.(P5|) still smaller). 

By a similar procedure, the projector for extracting the /c-resonance /^from the set L„ could 
be obtained: 

Qk{n) = e'^'^--f[(l + i^^). (24) 



The contamination of the v?fc(x) resonance could be more strong than for ip{x) due to the factor 
exp(/ifcr). 

It would be worthwhile to rewrite expression ( p3D as a form of Laplace representation: 

rco ^ ( p \ ~ 

if(x) = / e~P'^TT 1 ^(p,x)dp; where $(r, x) = / e~^'^^(p,x)dp (25) 

Jo fJi \ HiJ Jo 

and then to employ the identity 

n ( 1 - — I = exp ('-psi -p^^ - . . . -p"— - . . .V where Sfc = X] ( — ) • (26) 



We are grateful to A. V. Radyushkin who suggested us to use just the same form of the projector 
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Now, all the information about the spectrum is accumulated in Sj-parameters in the final 
formula: 

ip{x) = exp l^-pr - psi - JIp^J^ ^) + 0{e-^"+'^). (27) 



Expression (^71) is the main result of this subsection: it is very convenient for applications of 
any models of the spectrum {/ij}, and allows one to use different approximations. 

4.2 The spectrum model and the ground state WF 

We suggest here the following spectrum model: 

The positions of the first two resonances are /ii = 1.7, /i2 ~ 3 {GeV"^) in accord with exper- 
imental data and estimations of subsection 3.3. For subsequent /ij we adopt the "physically 
appealing" oscillator's model with a finite spectrum: fi^ = fi, /i4 = yU+l,... /i„ = /i + n — 3. 

Now we rewrite expression (^) in a form without the first two resonances, and all the model 
suggestions concentrate in the intergrand exponential: 

ip{x)= + (^1 + ^^ J\xp(^-pT-psi~Y.^p''j^^{p,T)dp + 0{^^ (28) 
where 

Our model allows one to apply any approximation to the integral in the r.h.s.. We shall not here 
describe the details (see Appendix A). It results in a "three-resonance formula" for expression 
(|28|) with a shifted argument r 

ifix) ^ $3^(3;, t)=(i + —] (1 + + -] $(r, x) = -. (29) 

V t^lj \ t^2j \ f^J Uf^i 

Note here that the shift r r + si effectively reduces the available region of the parameter 
M^: 

1 + SiM2 

Evidently, this "scaling" effect improves the stability of (p{x) with respect to M^. 

The "reduced" expression (p9|) depends on two parameters, fi and n. We can vary them 
to obtain the best region of stability for if{x) with respect to M^. The best fit leads to the 
values: ~ 4 GeV^, n = 11 for ansatz A2. For these values of parameters the form of the 
curve (f{x) doesn't depend on visibly in the region e (0.5 — 2) GeV^, and has a rather 
slow evolution when increases. The form of curve (in the stability region) is similar to the 
asymptotic WF again (see FigJ^), and the second moment of WF (^^) = 0.16 for = 1 GeV"^. 
So, we may conclude that the results of the "differential method" confirm the results for pion 
WF obtained by the "integration method" in Sec. ^ 
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5 Conclusion 



In this paper, we have considered a model sum rule for the pion wave function and pseudoscalar 
resonances based on the non-diagonal correlator The theoretical side (r.h.s.) of this sum 
rule depends only on the non-local condensate. We suggest two different ways to extract the 
wave functions of the ground and resonance states of the meson from the sum rule. Two kinds 
of ansatzes with different analytic properties for the condensate distribution have been used 
for final numerical estimates. We have obtained qualitatively similar results for the pion wave 
function by using both these methods and different condensate ansatzes: their forms are close 
to the asymptotic wave function. Moreover, we obtained the shape of the wave function and 
the position of the first resonance near mf ^ 1.8 GeV'^, {rn\^^p = 1.7 GeV"^), based on the 
first kind ( A2 ) of the ansatz. Our results confirm the qualitative conclusions about shapes of 
the wave function of pion and the resonance obtained in |jl|]. 
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Appendix 



A Three-resonance formula 

The oscillator spectrum model leads to the evident expression for coefficients Sk in (pBI): 



1 



To derive the three-resonance formula (p9|) we shall use a set of successive approximations ( [A. 2 
- A. 5). Substituting ( [A .11 ) into the identity ( P^D we obtain: 



k=2 k=2 



^ [/i ■ (1 - z,) ln(l - zi) - + 1) ■ (1 - Z2) ln(l - ^2)] |.,<i, (A.2) 
where Zi = p/ fi, z^ = p/(r;, + 1). 



exp[-(/i-(l-Zi)ln(l-;^i)-n-(l-Z2)ln(l-;^2))] ^ ^ ^^(^. (A.3) 



Simple numerical analysis demonstrates the validity of the approximate equation for the r.h.s. 
ot(H): 

"("-rt— 7^ (l-^)- (A.4) 



The main condition for the next approximation is that the Laplace representation r) for 
A2 is concentrated on a bounded support of about few GeV^ ^ and the integral in (^) over p is 
saturated in the interval (0,/i). So, we may hope that 



dp 



<l>(p, x) exp (-pr - psi) ^ f 1 + ^ j $(r, x) |r-.r+si • (A. 5) 



Substituting the set of approximations ( |A.2| - [A.5| ) into expression (|^), we arrive at three- 
resonance formula for '^{x)-,^ 

^{x\ ^ (^1 + £^ (^1 + £ j (^1 + ^ j $(r, a;) |,^,+,, . (A.6) 
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Figure 1: ^x,M^) with = 0.4 GeV\\me), = 0.6 Ge^^ (small-dashed), = 
1.0 Gey ^(big-dashed) vs. asymptotic pion wave function (dotted) 



Figure 2: A2-ansatz: ^n{u!,x) with N = l(line), N = 2 (small-dashed), N = 3(big-dashed) vs. 
asymptotic pion wave function (dotted) 



Figure 3: A34-ansatz: ^n{uj,x) with N — 6(hne), N — 7(small-dashed) vs. asymptotic pion 
wave function (dotted) 



Figure 4: A2-ansatz: $i(a;, 0.5) - ^2{uJ, 0.5) 



Figure 5: A2-ansatz: (pi{x,u;) with u; — 2.5 Gey^(line), u; — 3.0 Gey^(dashed) and u; — 
3.5 Ge\/2 (dotted) 



Figure 6: A2-ansatz: <^3r{x,j^) with = 0.5 GeV^{\me), = 1.0 Ge^^ (small-dashed), 
— 1.5 Ge\^^ (big-dashed) vs. asymptotic pion wave function (dotted) 
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